Real-time propagators at finite temperature and chemical potential 
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We derive a form of spectral representations for all bosonic and fermionic propagators in the 
real-time formulation of field theory at finite temperature and chemical potential. Besides being 
simple and symmetrical between the bosonic and the fermionic types, these representations depend 
explicitly on analytic functions only. This last property allows a simple evaluation of loop integrals 
in the energy variables over propagators in this form, even in presence of chemical potentials, which 
is not possible over their conventional form. 

PACS numbers: ll.10.Wx 



I. INTRODUCTION 

An advantage of the real-time formulation of field 
theory in medium over its imaginary-time counterpart 
is its close resemblance with the field theory in vac- 
uum P, 0, H, H|. A matrix propagator in this formu- 
lation is given, roughly speaking, by the vacuum propa- 
gator and another piece representing the medium. Also 
one deals here with an energy integral instead of a fre- 
quency sum. Thus, even though its matrix structure is 
looked upon as a source of complication, the real-time 
formulation has found many applications to study the 
properties of hadrons propagating through a medium 
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However, a difficulty arises with the conventional form 
of propagators in presence of one or more chemical po- 
tentials. Such propagators have theta-functions in the 
energy variable separating the contributions of particles 
and antiparticles in the medium. Then, unlike the case in 
vacuum, the energy integral in a loop contribution cannot 
be carried out by the simple contour method [HI, [l3[ ■ 

In this work we rewrite all the real-time propagators in 
terms of their Kallcn-Lchmann spectral representations 
in a suitable form. Besides the spectral densi- 
ties, these representations depend only on analytic func- 
tions, even in presence of chemical potentials, allowing 
contour evaluation of integrals in loop energies. The in- 
tegrations over the (known) spectral densities can then be 
performed, bringing out the different distribution func- 
tions for particles and antiparticles, when chemical po- 
tentials are present. 

As expected, our representations are also valid for com- 
plete propagators, in fact, for any two-point correlation 
functions. These are found to be only of two types, de- 
pending on the bosonic or the fermionic character of the 
operators in the two-point functions. The two types of 



matrix propagators are rather simple and enjoy a high 
symmetry between them. 

We begin with representations for the three- 
dimensional (spatial) Fourier transforms of the two-point 
functions. After selecting the time contour, we carry out 
the remaining Fourier transform with respect to time to 
get the desired representation. These, of course, repro- 
duce the conventional form of free propagators, once the 
integrations over their known spectral densities are car- 
ried out. 

In Sec 2 we consider the spectral representations for 
the spatial Fourier transform of the free propagators. In 
Sec 3 we show that the same form is also valid for the 
complete propagators. Then in Sec 4 we obtain the four 
dimensional transforms. In Sec 5 we work out an example 
of a loop integral with bosonic and fermionic propagators. 
Our discussion is contained in Sec 6. 



II. FREE THEORY 

Let us denote the ensemble average of an operator O 
by (O), 



(O) = Tr [e 



0]/Z, Z = Tr[e 



-/3(H-»Nh 



(1) 

where H and N denote the Hamiltonian and the number 
operator of the system under consideration, maintained 
at temperature 1//3 and chemical potential /x. 

We begin with the free propagators for a scalar field 
4>{x) (with no chemical potential) and a (Dirac) spinor 
field ip(x) (with chemical potential fx), 



D (o) ( 

S (0)( 



- x') = i{T c cf){x)4>{x')), 

— x') = i(T c ijj{x)ip(x')} 



(2) 
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where x M = (t,x). Here r and r 1 are any two 'times' 
on a contour in the complex time plane. The ensemble 
average with the Boltzmann weight and the consequent 
analyticity domain, —j3 < Im(r — r') < 0, constrains the 
contour to begin at — T, say, on the real axis and end 



2 



at — T — if), nowhere moving upwards [16|. Apart from 
these requirements, we let the contour be arbitrary at 
this stage. Here T c (and 6 C below) denote time ordering 
(and 9- function) with respect to such a contour. 

Our form for the propagators in momentum space fol- 
low directly from the following spectral representations 
of their spatial Fourier transforms [HI , 



+ DC 



dk' 



2tt 



V° ) (*k*)e- tt&(T - r ' ) {0c(7 



/(fto)K3) 



5(°)(r-r',p) 



dp'p 
2ir 



a^(p' o ,p)e-^-^{0 c (r - r') - /(p{,)}(4) 



where p^ and are the bosonic and the fermionic 
spectral functions for the free theory, 

p m (k ,k) = 27re(k )S(k 2 -m%) (5) 
aW(po,p) = 2ire( Po )(]b + m F )5( P 2 - m%) (6) 



and the functions / and / are given by 
1 



/(fco) = 



,/3fco — 1 



/(PO) 



1 



e /3(p -M) _|- 1 



As fco and po run over the entire real axis, / and / cannot 
be interpreted as particle distribution functions, which 
are 



"(fco) 
n+(po) 
n-(po) 



e /3|fcol — 1 
1 



e /3(M-p) -f 1 
1 



(boson) 

(fcrmion) 



(antifcrmion) (8) 



e P(.\Po\+v) -|- i 

But they may be readily expressed in terms of the latter 



ones as 



/(fco) = ne(k ) - 6(-k ), f( Po ) = JV l£ (po) + 0(-po) 

(9) 

where 



Ni( Po ) 



APo) 



A-Po) 



(10) 



A simple way to derive the results (3) and (4) is to 
note that each of the propagators in medium satisfy the 
same differential equation as the corresponding one in 
vacuum. Taking their spatial Fourier transforms we get 
one dimcntional equations in time r. A particular so- 
lution to such an inhomogeneous equation is given by 
the one with Feynman boundary condition. To get the 
most general solution, we add to it the two independent 
(plane wave) solutions of the homogeneous equation with 
arbitrary coefficients. When these coefficients are deter- 
mined by applying the thermal, so-called Kubo-Martin- 
Schwinger (KMS) boundary condition [UGl], we § e t the 
above solutions. One can, of course, readily verify these 
solutions by showing that they satisfy the respective dif- 
ferential equation as well the boundary condtion. 



III. INTERACTING THEORY 

Before we proceed, let us consider complete propaga- 
tors for interacting theories in the Heisenberg picture and 
obtain the Kallen-Lehmann representation for their spa- 
tial Fourier transforms. Consider first the full propagator 
for the interacting scalar field $(x), 



D(x - x') = i(T c $(x)$(x')) 



(11) 



Now it does not satisfy a simple differential equation any 
more. But we may analyze the field product by intro- 
ducing a complete set of eigenstates \m) of the four- 
momentum operator with eigenvalues {E m ,k m ). We 
insert this set twice in Eq. (11), once to evaluate the en- 
semble trace and then to separate the operators in the 
product. Using translational invariance to extract the 
x-dependence of the operators and introducing a delta- 
function in the energy variables, the spatial Fourier trans- 
form of (&(x)Q(x')) in the T c -product becomes 



cP xe - i *<*- 3 ")(&(x)&(rf))= I ^-e- tk '^ T - T,) M + (k' ,k) 



(7) where the spectral density M + is given by 



(12) 



M + (k ,k) = (2tt) 4 ^ e - 0Em S^(k + k rn -k n )x 

(m|*(0)|n)(n|$(0)|m)/Z(13) 

In just the same way we work out the Fourier transform 
(§(x')&(x)) with a spectral density M~, that may be 
obtained from M + by changing the signs of fc m and k n in 
the delta-function. The two spectral densities are related 
by the KMS condition 



M+{k) = e fjk< >M-(k) 



(14) 



which may be obtained simply by interchanging the 
dummy indices m, n in one of M ± (k) and using the en- 
ergy conserving 5-function [l9l ]. 

We now define a third spectral density related to the 
commutator of the field operators, 

p(k) = M+(k)-M-(k) = J d A xe lk < x - x '\[^{x)^(x')]) 

(15) 

in terms of which (k) can be expressed as 

M+(k) = {1 + f(k )}p(k), M-(k) = f(k )p(k) (16) 

Using these relations, we see that the three-dimensional 
Fourier transform of the complete propagator is given 
exactly by Eq.(3) with the free density replaced by 
the interacting one, p. 

The complete propagator for the interacting spinor 
field #0), 



S(x-x') = i{T c V(x)V(x')} 



(17) 
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can be treated in the same way. The three dimensional 
Fourier transform of (^(x) 1 $(x 1 )) is given by an equation 
like (12) with the spectral density 

M + (p ,p) = (2^) 4 J2 e-^-^5^ip + p m -p n ) x 

(m|*(0)|n)(n|*(0)|m)/Z, (18) 

while that of (^f(x')^f(x)) is given by a second spectral 
density M~ (po-,p), obtained from the first one by inter- 
changing ^(0) and ^(0) and changing the signs of p m 
and p n . Here the integers N m are the eigenvalues of the 
operator N, the fermion number conservation restricting 
them to N m = N n — 1. The two spectral densities are 
again related by the KMS condition 

M+(p) = e^°-^M-{p) , (19) 

Defining a third spectral density in terms of the anti- 
commutator of the field operators, 

a{p) = M+(p)+M-(p) ee J d 4 xe lp - {x - x ">({^(x), *(x')}) 

(20) 

we can express the first two as 
M+(p) = {l-f(p )}a(p), M"(p) - f( Po )a(p) (21) 

We thus again get the complete spinor propagator as 
given by Eq.(4) with er ) replaced by a. 

It should be noted that in deriving these representa- 
tions for the complete scalar and spinor propagators, we 
have not used any property specific to these field opera- 
tors, except to define their spectral densities by the com- 
mutator and the anticommutator of the respective fields. 
We thus conclude that these are the most general repre- 
sentations for the spatial Fourier transforms of two-point 
functions of any bosonic and fermionic local operators. 
Below we construct their real time representations, again 
without requiring any property of the spectral densities 
at all, so that although we refer to them as propagators, 
they actually remain valid for any two-point function. 

In passing, let us compare the derivation of the Kallcn- 
Lchmann representation for the ensemble average of two 
time-ordered fields with that for its vacuum expectation 
value. In the vacuum case, the intermediate energies arc 
all positive and Lorentz invariance makes the spectral 
densities to depend on 6(ko) and fc 2 (= k$ — k 2 ) (6(p ) 
and p 2 {= Pq — p 2 )). There is, of course, no KMS con- 
dition here, but since the commutator (anticommutator) 
vanishes outside the light-cone, the two parts of spectral 
densities must be equal [20]. In the medium, there is no 
Lorentz invariance (or we may have it at the cost of in- 
troducing a four- velocity vector for the medium). As a 
result, the vanishing of the commutator (anticommuta- 
tor) does not lead to any simple condition on the spectral 
densities. 
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FIG. 1: Contour in time plane for real time formalism 



IV. FORM OF PROPAGATORS 

Having established the results for the spatial Fourier 
transforms of the propagators, where we kept their time 
coordinates on an arbitrary contour in the complex time 
plane, we now choose an appropriate contour to get the 
real time field theory. Our choice is that of Fig. 1, which 
for T — > oo, reduces to two parallel lines, the real line and 
the one shifted by —1(3/2, to be denoted by subscripts 1 
and 2 respectively [l[. Thus n = t, r 2 = t — i/3/2, 8 c (ti — 
t[) = 9(t—t'), c (ti — t 2 ) = 0, etc. We can now define the 
Fourier transform of the four components of a propagator 
with respect to real time. We thus get from Eqs. (3, 4) the 
four dimensional Fourier transform of the components 
respectively of the boson and the fermion propagator as 

/+00 
dte iko ^-^D(T a -Tl,k), 0,6 = 1,2 
-OO 

(22) 

/+00 
dte^-^SiTa-T^p), a, b =1,2 
-OO 

(23) 

It is now simple to work out the matrix propagators. 
For the scalar field, it is 

f°° dV 

D(k ,k)=j ^p(k' ,k)A(k' ,k ) (24) 
with the elements of the matrix A given by 

A n = -A* 22 = \ + 2mf(k' )S(k' Q - ko) 

Kq Kq 11] 

Aia = A21 = 27ne^ 2 f(k' )S(k / a - fc ) (25) 
The corresponding results for the spinor propagator is 

S( Po> p) = f°° ^a(p' ,m(p' ,p ) (26) 

where the matrix $7 has the elements, 

n n = -n* 2 = : 2nif(p' )S(p' -po) 

12 = -e^fiai = -27ne/Po/ 2 f(p' Q )S(p' - p ) (27) 

Eqs. (24,25) and (26,27) constitute our form for the spec- 
tral representation of any bosonic and any fermionic two- 
point function respectively. The two matrices A and 
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$1 are symmetrical, being related by the replacement 
/ <-» — /. These describe the common propagation prop- 
erties of all elementary (particle) and composite fields, 
the dependence on spin, etc being contained in the re- 
spective spectral function. 

We also note that the distribution-like functions / and 
/ given by Eq. (7) which occur in Eq.(25) and (27) are an- 
alytic functions, in contrast to the distribution functions 
themselves given by Eq.(8). In particular, the matrices 
A and do not contain any (9-function in the energy 
variables fco(po), which are present in the conventional 
form of the propagators in presence of chemical poten- 
tial. Thus, as we show in Sec. 5, the use of the above 
spectral reprsentations avoids the difficulty encountered 
in integration of loop energies over the conventional form 
of propagators. 

Recalling Eq. (9), we can, of course, readily restate the 
matrix elements given by Eqs. (25, 27) in terms of the 
particle distribution functions themselves as 



An 



_jY* = 

22 k' - k ~ iT]e(k Q ) 



2irine(ko)5(k — ko) 



22 Po -Po ~ w?e(po) 
fi u = -e^fi 21 = -2TrieW 2 N 2 e(p )S(p' a - p ) (28) 

where Ni is defined by Eq. (10) above and N 2 by 



N 2 (po) = y/n+{l-n + )0(po) - y/n-{l - n-)0(-po) 

(29) 

Note that we choose primed variables (k' ,p' ) for the 
functions multiplying the delta-functions in Eqs. (25,27), 
while we choose unprimed ones (k ,p ) for such func- 
tions in Eq. (28). The first choice allows us to integrate 
over loop energies in a simple way (see Eq. (46) below). 
The second choice makes the matrices diagonalising the 
propagators independent of the integration variable in 
the spectral representations (see Eqs. (32,33) below). 

The matrices A and as defined by Eq. (28) can now 
be diagonaliscd respectively by U and V , 



U(k Q ) = 
V( Po ) = 



+ n v 
n 



(30) 



-0fi/2 AT-/./XT I 



which in turn diagonalises the two-point correlation func- 
tions themselves, 



D(k ,k) = U(k 



D 



JL* ) U(k ) (32) 



1 o -T ■ 



where the analytic amplitudes D and S are given by 



D(k ,k) 



dk' n 



p(K,k) 



2ir k' — ko — ii]e(ko) 



(34) 



S(po,P) 



d Po vjp'oiP) 

2tt Po - po - ^e(Po) 



(35) 



To find the spectral densities p and a, we need not 
refer to their defining equations (15) and (20) any more. 
It suffices to calculate the imaginary part of any one, say 
the 11-component of the two point functions to get them. 
Thus from Eqs. (25) and (27), we get 

p(k ,k) = 2tanh(/3/co/2)ImL>n(fco,fc) (36) 
^(po,^ = 2coth(/3(p - M )/2)Im5n(p ,p) (37) 

We close our discussion of the form of two-point func- 
tions by obtaining explicitly the conventional form of the 
free propagators for the scalar and the spinor fields @, 0] . 
With the free field spectral densities given by Eqs. (5, 6), 
Eqs. (34, 35) give immediately 



"1 



= A(k,m B ) 



S ° (p) = {i> + m F) A(p, m F ) 



(38) 



A12 = A 2 i = 2-Ki\Jn{\ + n)e(k )5(k' - ko) 
fin = -O22 = — — — ; — 7 - 27ri7V 1 e(p )(5(p -p ) 



which are just the Feynman propagators in vacuum. 
Then multiplying out the matrices in Eqs. (32, 33) we 
get the scalar propagator as 



D (0) (k ,k) 



dn d 12 
d 2 \ d 22 



(39) 



where the matrix elements are 



(in = — d 22 = A(fc, tub) + 2ninS(k 2 — m B ) 

dn = d 21 = 2iriy / n{l + n)6(k 2 - m|) (40) 



and the spinor propagator as 



S^°'(po,p) = + m F ) 



with the matrix elements 



Sll S12 

s 2 i s 22 



(41) 



Sll 

S12 



s 22 = A(p, mp) — 2niNiS(p 2 — trip) 

-2Trie^ /2 N 2 5{p 2 - m 2 F ) (42) 



-e^32i 



V. AN EXAMPLE 

Finally we evaluate in some detail the energy integral 
in a loop involving a scalar (with p, = 0) and a spinor 
(with n 7^ 0) propagator, both exact, using their spectral 
representations given by Eqs. (24-27). The matrix ampli- 
tude of such a loop, call it S, belongs to the complete 
spinor two-point function and by Eqs. (35,37) its analytic 
amplitude S must have the form 



£(po,p) 



d Po <n(Po 1 P) 
2ir p' -p - iv^Po) 



(43) 
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where the one-loop spectral density is given by Eq. (37) 



a.s 



ai(p ,p) = 2coth{/3(p - M)/2}ImEii(po,£). (44) 

Ignoring the spin and isospin structures of the interac- 
tion vertices, the 11-component of the matrix amplitude 
is given by 



p(qZ,f- q){l f(q' ) + /(<#)}*(?„ -q' - <$) 

(49) 

Further, if the propagators are the free ones, one can 
easily work out the q and g ' integrals over the mass 
shell (5-functions in the spectral densities to get [ cJi = 
yjq 2 + rri 2 F , uj 2 = \/(p-q) 2 +m A B } 



/d 4 q f d?q 



f^Ul P ^p(q'^-q)I(q' ^,Po) (45) 



where 



f+oo ^ 

I (q'o,q'o,Pa) = i I ■^jr^ii(?o>9o)An(g'o,Po - Qo) 



Rewriting the propagators in / as 

i - f(q'o) 



(46) 



fWo) 



ftn(q' ,qo) = — 

%-Qo- "7 Qo - 9o + iV 

An(M = 2±ML.^M^ (47) 

%-Qo- irj q -qo + iv 

we can integrate the resulting four terms over go by just 
closing the contour in the go plane to get 



I{q'o,q'o,Po) 



{l-/(?o)Hl + /(So)} , /(3o)/(So) 



Pq - q'o - Qo + ir \ Po - <2o - Qo - ir \ 

(48) 

Had we taken the spinor propagator from Eq.(28), the 
^-functions in it would not have allowed such a simple 
evaluation of 7 [T^. 

The imaginary part of En is given by that of / 

Im/(g ,g ',po) = -n{(l-f)(l+ f)-ff)}S(j> -q -q^) , 

which has the 'wrong' relative sign between the two terms 
in the bracket. It can be reversed by taking out a factor 
of tanh{/3(po — m)/2}- Then F reduces to 



0i(po,p) 



-2tt 



q 



c(?o><7) 



(2tt) 3 2tt ~ ^"'^ ./_™ 2tt 



dq'o 



(2tt) 3 

{(1 - ?i + + n)5(p - uj 1 - lu 2 ) + (n+ + n)<5(p - &i + w 2 )} 
— (wi — > — wi, W2 — ► — W2,n + — ► n_) ] , (50) 

in agreement with the computation in the imaginary time 
formulation plj . 

VI. DISCUSSION 

Despite many applications of the real-time propaga- 
tors in medium, their conventional matrix form is not 
satisfactory in that the matrix components contain non- 
analytic functions, particularly when chemical potential 
is present. This causes difficulty in carrying out integrals 
in loop energies over such form of propagators. 

In looking for improvement on this point, we take our 
cue from the Kallen-Lehmann representation for the vac- 
uum two-point functions, where the propagation proper- 
ties common to all fields are displayed explicitly by ana- 
lytic functions, keeping dependence on spin etc. confined 
to the spectral densities. We are thus led to study this 
representation at finite temperature and chemical poten- 
tial for all two-point correlation functions. 

Indeed, we find spectral representations here to de- 
pend explicitly only on analytic, distribution-like func- 
tions, along with propagation function expected in vac- 
uum. They arc only of two types, one for the bosonic 
and another for fermionic two-point functions. Also 
they are simple and symmetrical between the two types. 
We choose this representation even for free propagators. 
Then the loop integrals in the energy variables can be 
easily carried out, irrespective of the presence of chemical 
potentials in them. Integrals over the spectral densities 
can then be performed. 
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